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The problem of optimally estimating an unknown unitary quantum operation with the aid of 
entanglement is addressed. The idea is to prepare an entangled pair, apply the unknown unitary 
to one of the two parts, and then measure the joint output state. This measurement could be an 
entangled one or it could be separable (e.g., measurements which can be implemented with local 
operations and classical comunication or LOCC). A comparison is made between these possibilities 
and it is shown that by using nonseparable measurements one can improve the accuracy of the 
estimation by a factor of 2(d + 1) jd where d is the dimension of the Hilbert space on which U acts. 



PACS numbers: 03.67.-a 



I. INTRODUCTION 

Consider a one-qubit unitary gate, the following ques- 
tion arises: "how to characterize it?" This question is 
motivated by recent experiments in quantum optics pj. 
A possible approach is to prepare many known states and 
use them as inputs, and then measure the outputs that 
they produce; this is known as quantum process tomogra- 
phy |2j. It turns out that one needs as inputs a basis of 
the Hilbert space plus some linear combinations thereof. 
The disadvantage of this approach is that, in many prac- 
tical situations, such a set of states is not feasible in the 
laboratory pj. 

Another strategy is described in Refs. 0, 0, Q It is 
enough to use a single bipartite entangled state; one of 
the states is used as input for the quantum operation and 
nothing is done to the other one, then the two qubits are 
measured, as shown in Fig. ^ 
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FIG. 1: The use of a single entangled input state suffices. 

In Ref. it is pointed out that in this setup there is a 
one to one correspondence between the quantum opera- 
tion and the joint output state. A maximally entangled 
state is used as input and then the three components of 
the spin in both output particles are measured. One can 
ask whether it is possible to find a more accurate mea- 
surement. Also, is it possible to find a measurement that 
performs as well as the one in Ref. which has less out- 
comes? It will turn out that the answer is that one can 



find a more accurate measurement but this measurement 
is nonseparable. It is also possible to find a measurement 
with less outcomes. 

In Ref. [U it is proven that a maximally entangled 
pure state is a good input state in the sense that if |</>) £ 
C 2 <£> C 2 is maximally entangled, then 

VponC 2 ®C 2 , H p {0) < H WW (J9) 

where H p (9q) is the quantum Fisher information matrix 
(QFI) for the joint output state, at 9 — 9q, if the input 
state is p. This quantity is defined for example in Ref. 
and explained in more detail in the following section. 
The inverse of this matrix is a lower bound (quantum 
Cramer-Rao bound or QCRB) for the mean-square er- 
ror of estimators based on arbitrary measurements of the 
output state. A maximally entangled state is a good in- 
put also in the sense that the QCRB can be achieved if 
and only if the input state is maximally entangled. The 
problem here is that, as will be shown later, the mea- 
surement that achieves the bound is actually a basis of 
projectors onto maximally entangled states. This mea- 
surement can be performed using nonlinear optics but is 
far from being standard. 

But perhaps the improvement in the estimation 
through the use of entangled measurements is not very 
large. Are entangled measurements worth the trouble? 
The main aim of this paper is to show that the answer 
to this question is positive. The value of entangled mea- 
surements will be quantified precisely. 

Before continuing with this discussion it is necessary to 
explain what is meant with "better" and "best" and how 
the quality of different positive operator valued measures 
or POVMs is actually going to be quantified. 



II. QFI, FISHER INFORMATION, AND QCRB 
A. Quantum Fisher information 
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Suppose that the quantum state density matrix a on 
C d is parametrized by 9 £ O C R p where p is the number 
of parameters (less than or equal to d 2 — 1 for mixed 
states, 2d — 2 for pure states). In our case a would be 
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the joint output state. Define the symmetric logarithmic 
derivatives Ai , . . . , \ p as the self-adjoint operators that 
satisfy 

<r,i(j9) = d ei a{6) = \[a(B)\{6) + \i(9)a(e)\. 

For pure states, a = \ip)(ip\^ they simply are Ai = 2(7,,. 
The QFI is defined as the p x p matrix with elements 

HijiO) =Rc tr[a(0)A l (0)A,(0)] 

which for pure states reduces to 

H ij (e)=Re(k(9)\lM) 
where \k(6)) ^ \ (9)\^(0)) . 

B. (Classical) Fisher information 

Take a POVM with elements Mi, ... , M n . The Fisher 
information matrix (FI) for this measurement is the pxp 
matrix with elements 



tv[p(e)M £ ] 



For an estimator 9 and a measurement M, locally unbi- 
ased at 9o 1 , the (classical) Cramer- Rao bound is satisfied 

v(M,e 0> 9)>i(M,e )- 1 , 

i.e., the FI is the smallest variance that a locally unbiased 
estimator based on this measurement can have. This also 
means that if one of the eigenvalues of I is zero, then the 
variance of the function of the parameters correspond- 
ing to that eigenvalue is infinity and therefore cannot be 
estimated. 

If one has N copies of the quantum state and performs 
the same measurement on each of the copies then the 
FI of the N copies, I N , satisfies I N (M,9) = NI(M,9) 
where I(M, 9) is the FI of one system. It follows that 

V N (M, 9q, 9) > I N (M, 9 )~ 1 = I(M,e )- l /N. 

It is a well known fact in mathematical statistics that 
the maximum likelihood estimator (MLE) in the limit 
of large N is asymptotically unbiased and saturates the 
classical Cramer-Rao bound. Moreover no other reason- 
able estimator (unbiased or not) can do better. 



C. QCRB 

The QCRB states that for any measurement M 

I(M,9)<H(9). (1) 

In other words, H(9) — I(M, 9) is a positive semidefinite 
matrix. 

This bound is not achievable in general. A theorem due 
to Matsumoto states that for pure states, the bound 
is achievable at 9 = 9q if and only if 



Im(*i(0o)M0o)> =0. 



(2) 



Furthermore, if condition @ holds, there is a measure- 
ment with p + 2 elements that achieves the bound. 

In analogy with Q, measurements will be compared 
using the quantity 

tr H{6)- l I{M, 9) 

which is always less than or equal to p, the number of 
parameters. For example, for the measurement used in 
111, trH- 1 ! = 1. 

One needs to use a quantity like this because of the 
extra complexity that quantum theory adds to the prob- 
lem. Namely, in the most general case there is no 
POVM that achieves equality in QJ. Typically, for any 
POVM M\ which cannot be improved, one can find an- 
other POVM M 2 such that neither I (Mi, 6) < I(M 2 ,9) 
nor J(Mi,0) > I(M 2 ,6) are satisfied. The bound Q 
is sharp, i.e., H(9) is the smallest matrix larger than 
I(M, 9) for every M. The difficulties can be overcome by 
using a single number (instead of a matrix) to quantify 
the performance of a POVM. This defines an achievable 
bound and any two POVMs can be compared according 
to this quantity. Of course, no single number can be an 
absolute quantification of the performance of a POVM. In 
applications one must decide what one wants to estimate 
and accordingly assign weights to the mean square error 
of the parameters to be estimated. This comes down to 
using a quantity such as trG{6)I{M,6). One needs to 
maximize this quantity for a general 2 G{9) > tailored 
to one's specific needs. In this paper the general problem 
is not solved. Only the case G(9) — H(9)~ 1 is considered. 
There are several good reasons for this choice: 

1. Since H(9) is the smallest upper bound for all the 
I(M, 9), it defines a natural scale in which to com- 
pare them. 

2. tr H(9)^ 1 1(M, 9) is parametrization invariant. 



2 This quantity still has the property that if the inequality 

1 This means that the expectation of the estimator satisfies tr G(6»)7(Mi, 0) > tr G(9)I(M 2 , 9) holds, then I(Mi,$) jt 

E M ,e o 0i) = 9 0i and dg E M ,$0i) I = kj- I(M 2 ,9). 

I 6— 8 
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3. H (9) is closely related to the fidelity between true 
and estimated output states: the metric generated 
by H (9) is locally identical (up to a factor of 4) to 
the Bures distance, enures (p, c) 2 = 2(1— ^JTifi, a)), 
where T is he fidelity which for pure states can be 
defined as T(\^\MW\) = I(V#)I 2 - 

4. Finally, the use of this quantity allows one to obtain 
simple and striking results. 



III. THE CASE d = 2 

A. Entangled measurements 

In H it was shown that in this case the best input is 
any maximally entangled state. Here the singlet state 
\t) = [| 10) - |01)]/\/2 will be used. The output is 
then \i/}(a, 9, 4>)) = [U(a, 6, 4>] <g> l)|r) where U(a, 0, <f)) = 
cos a 1 + i sin a ng^-a is a 2 x 2 unitary matrix, ftg^ is the 
unit vector (sin#cos</),sin(9sin</>, cos (9) parametrized by 
its polar coordinates, and a = (ci, 02, 03) are the Pauli 
matrices. 

It is quite straightforward to calculate that the QFI is 

/ 1 

H(a,6,(f)) = 4 (sin a) 2 

\0 (sin a sin6») 2 

and as expected lm(li(a,8,<p)\lj(a,6,4>)) = 0. One can 
find a simple measurement that achieves the bound; in 
fact, any measurement of the type 

M a = \b a )(b a \ a=\,...,p+l, 

m+l 

M p+2 = t - J2 M a, 

£J (3) 
\b a ) = 2_,o af3 \mp), 

(3=1 

\m k ) =J2( H ~ 1/2 hi\k), K +1 ) = |0), 

with oa (p+l)x + real orthogonal matrix satisfying 
o a ,p+i 7^ achieves the bound. For example measuring 
the Bell basis 



Bell [00> - <00| — <11| 
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f Bell_ |00) + |11) (00| + (11| 

M BeU = 



V2 V2 ' 

M |01) + |10) (01| + (10| 



(4) 
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achieves I(M Belt ,6) = H{9) for all 6 and therefore satis- 
fies 



everywhere; this is three times the value achieved in Q). 
This measurement, which has been implemented using 
non-linear optics Q, is not widely available in quantum 
optics labs. On the other hand, a POVM with the two 
components 



Mi 
M 2 



M fe Be ", 
1 - M, 



Bell 



(6) 



for k = 1, 2, 3 or 4 has been implemented with linear op- 
tics and is much more standard than measuring the whole 
basis. This is a POVM with only two outcomes (which 
might be an advantage for its practical implementation) 
and calculation shows that it satisfies tri? -1 / = 1 ev- 
erywhere, which is as good as the measurement in [J. 
It will be shown in the next section that it is actually 
as good as any separable measurement (in terms of the 
value of tr H~ l I). This measurement does have a serious 
drawback, namely that one can only identify one func- 
tion of a, j3 and 4>. The drawback can be overcome, for 
example, by measuring © for fc = 1, 2 and 3 each in one 
third of the available copies. In this way one is able to 
identify all three parameters, tri? -1 / = 1 is still satis- 
fied, and finally this new POVM should not be harder to 
implement than the previous one. Again here everything 
depends on what one wants to estimate. 
A measurement with three elements 



Mi 



M 2 = M; 



Bell 



(7) 
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1 - M, Be " - M 



Bell 



tr H~ L (6)1 (M Bell ,0) = 3 



(5) 



for some fc 7^ I has also been implemented with linear 
optics. In fact, it has been shown Q that, with linear 
optics, this is the best one can do. This POVM sat- 
isfies trff -1 / = 2. This is twice the value that can be 
achieved with any separable measurement. Note that this 
measurement has the same weakness as the previous one: 
it cannot identify all three parameters (it identifies two 
functions of a, 9 and 4>). It is easily possible to overcome 
this difficulty in a similar way as before. 



B. LOCC measurements 

How well can one estimate U using only LOCC meas- 
urements: measurements that can be implemented lo- 
cally and with the aid of classical communication be- 
tween the two parties. In fact, to begin with, the larger 
class of separable measurements will be studied: mea- 
surements whose elements are positive combinations of 
products of one dimensional projectors. This definition 
of "separable" is a slight generalization of that of [Toj . 
where it was shown that there exist separable measure- 
ments which are not LOCC. Nonseparable measurements 
are called entangled. 

Consider a separable POVM with elements = 
Ei^i(IVfi) ® ® (VfiD- One can refine this 
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POVM to obtain another POVM with elements that 
are proportional to one-dimensional projectors = 
<*W> ® ® (V-fJ). By relabeling U - £ 

one obtains 

M e = ce(|^)®|^f))((^|®(Vfl)- (8) 

The Fisher information corresponding to this refined 
POVM is greater than or equal to the Fisher information 
of the original POVM. Thus, since one wants to maximize 
the FI, one may restrict oneself to measurements of the 
type described in Eq. JSJ). 

For the calculations that follow it is more convenient 
to express Eq. JHJ using the Pauli matrices 



1 + 0£ • a t + b s -a 



with \a^\ = \b$ \ = 1 and > 0. The condition Y^t A^f 
1 can be rewritten as follows: 



6 e 

^2 c^b s =0,^2 c e a£kbti = 0. 
i 6 



(9) 



1 ® 1 - ^ °% ® a i J 



Since 

|10)-|01) <10| <01| _ 1 
y/2 v/2 4 

the density matrix of the system can be written as 

1 / 3 

p(a, 0,(f>) = - I 22 U{a, 6, ^a^U^a, 6, (f>) ® ct, 

V i=l 

The probabilities are then 



PS = 



C S 



a sr 



and i tr (UdiWaj) can be calculated to be 

3 

cos 2a <5y — sin 2a £ijk n k + 2 sin 2 a riiUj . 



fc=i 

Substituting this into the expression for the probabilities 
one gets 

= — ( 1 — cos 2a (a^ • 6g) + sin 2a (&£ x ■ n) 

■ V ^ ~ ~\ (10) 

—2 sin a (n • )(n • 6^) ) . 
After some not very interesting manipulations one finds 
1 ({Pi. a ? , (p^) 2 (^) 2 



\ 4 4 sin 2 a 4 sin 2 a sin 2 i 
= ^ ^1 + cos 2a (og • &j) — sin 2a (6^ x ■ n) (H) 

+2 sin 2 a (n • <%)(n ■ b^) 



Finally, using the conditions © one obtains that for sep- 
arable measurements of the type (JHJ) 



tr[H~ 1 (9)I(M, 9)] 



(j>i,4>? 



Pi V 



4 sin 2 a 4 sin 2 a sin 2 i 



(12) 



= 1. 



Any separable measurement can be refined to a mea- 
surement of the type JHJ. Therefore for all separable 
measurements M scp 

ixH(0)- 1 I(M„ Bp ,0) < 1. 

This bound therefore also holds for LOCC measurements. 
Since there are LOCC measurements of the type (JHJ), the 
bound is achievable with LOCC measurements. 



IV. THE CASE d > 2 
A. Entangled Measurements 

Before starting with any calculations it will be shown 
that the quantity that is being analyzed 

f(6) = sup tr H~ l (6)1(6, M) 

M 

docs not depend on 9. 

For any 9\ and 6*o there exists a unitary matrix V such 
that VU(9o) = U(9i). It is easy to see that for such a 
choice 

trH- 1 (6 1 )I(6 1 ,(V®l)M(V®±) f )=trH- 1 (6o)I(6 ,M) 
This implies that 

suptri/~ 1 (6»i)/(6»i,A/i) > snptTH- 1 (9 )I{9o,M ). 

Mi Mo 

Thus f(9\) > f(9o), but since 9q and 9\ are arbitrary, the 
function / must be constant. Therefore one can choose 
any value of the parameter to perform the calculations. 
One of the implications this has is that if one proves 
that the QCRB can (not) be achieved at one value of 
the parameter, then it can (not) be achieved everywhere 
(anywhere). 

In [j| it is mentioned that in dimension d > 2, it is 
no longer true that a maximally entangled state maxi- 
mizes the QFI; however it is still true that the QCRB is 
achieved if and only if the input state is maximally en- 
tangled. In order to prove the first statement it is enough 
to find a counter example. This is not difficult to do for 
example in d — 3. The second statement is also not dif- 
ficult to prove and because of the last discussion it will 
be enough to do it for U equal to the identity. 

A SU (d) matrix can be written as cxp ^i Y^L=i 6 a To^j ■ 
Here 9 € R d _1 and the T"s arc in the su(d) Lie Algebra. 
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They are traceless self-adjoint matrices and are chosen so 
that they also satisfy: 

tr(T a Tp) = S a0 . 

For U close to the identity (or 9 close to zero), 

d 2 ~i 



U Ps 1 + i e *T a 



a=l 

the input state can be written as Ylkl Rkl\kl)- Normal- 
ization implies tr RR^ = 1 where R is the d x d matrix 
with elements Rki- Since RR> has trace one and is self- 
adjoint it can be written RR^ = t/d + Yl a ^a.T a where 
the i's are real numbers. At the identity the output state 
satisfies 



\i>)=Y,Rki\kl), 
kl 

\il>, a ) = i^2RuT a \k)<g)\l)\ 



(13) 



kl 



the \l a ) vectors defined in section [Q] can be written as 

\l a ) = 2(\ta) + (t^M) 

and the condition for achieving the QCRB J2J) becomes 

2ti(RR*[T a ,T Y 



lm(l a \l0) = 4Im(V',a|^',/3) 



0(14) 



for all a and /?. This implies RR^ = l/e? 3 . This means 
that the input state is maximally entangled 4 . 

For the calculations the maximally entangled state, 
2fc=i \kk}/Vd, is used. H can be very easily calculated 
to be 

4 

H a (3 = ^af3- (15) 

Since the QCRB can be achieved 

sup tr H- 1 {8)1 (M, 9) = d 2 - 1. 

M 



B. LOCC measurements 

It will be shown here that for all separable measure- 
ments Ms e p the following holds 



tiH- 1 (9)I(M scp ,9) < 



d(d-l) 



(16) 



3 Since (RR^ — € su(n) and su(n) is a perfect Lie algebra 
(i.e. can be spanned by commutators), Eq. 1141 may be rewritten 

as vy esu(n ) tr[(flflt _ t/d)Y] = 0, this implies RRt - i = 
because the trace form is non-degenerate. 

4 The condition for a bipartite state to be maximally entangled is 
that the partial trace should be proportional to the identity. In 

tr 2 [ip)(i>\ = Rltf. 



This shows that if one allows nonseparable measure- 
ments, the estimation can be improved by a factor of 
2(d+ l)/d with respect to separable measurements. This 
is always more than twice. 

In order to prove Ea. (|16f) a particular representation 
for the T's will be chosen, namely: 



. s \k)(i\ + (-ir\i)(k\ 

T k is=i -j= k > I, s = {0,1}, 

d 

Tjn = ^ ' C m fc |fc) (fc| 



(17) 



m = 1, . . . , d — 1 , 



fe=i 



where the coefficients c m k obey 



^ Cmk = 0, 
fe=i 

d 

^ CmkCnk ^mn- 



(18) 



k=l 

From these two one can derive the relation 

d-l 

/] CmkCml = 5ki — -. 
m—1 

Measurements of the form 

M $ = cel^X^I = ce|o { )(o{| ® \h)(b s \ 

arc considered. The quantity of interest is 



trff" 1 / = - tr/ 
4 

= d v ((^|V», a ) (V>|^) + (^|^)(y>, a |^)) 2 



(19) 



So 



E 



c s 



Re 



<V#e> 



l(<^)l 5 
d 2 -i 



(20) 



a=l 



The second term in the previous equation is easy to cal- 
culate, 



d d2_1 d^^ 1 



(21) 



~ 2 2-< d 2 ' 

a=l 

but for the first term a little more work will be needed. 
One needs to calculate 



(hit 



our case 
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For a — {kls} 



,•8+1 



kls 



[<a#}(& e |0 + (-l)> € |0<«L 



E<^ 



(22) 



■( O€ |fc)(6 £ |0( O £|/)(6 € |fc). 



s=0 



Since the last expression is symmetric with respect to 
exchanging k with I, J2 k >i = \ Efc^ = j(Eh ~ Efc=i) 
and 

fe>; s=o fe=i 
In the case a = m 
d-i 

£<<^,™> 2 



m—1 



- <2 d-1 

-j £ £ c mk c mJ< a #>< 6 #)( a d i }( & dO 



fcj=l m=l 



d ^d 

kl 



Finally, substituting the previous equation and Eq. 1211) 
in (12 UH one obtains the desired result, namely, for any 
separable measurement M of the type © 



ixH -\e)i{M,e) = d -^—^ 



(25) 



Of course this impliesEq. I|16|) . The argument for LOCC 
measurements is the same as for the two dimensional case 
and one obtains the same bound for them. 



k=l 



V. CONCLUSIONS AND OPEN PROBLEMS 



In this paper it has been shown that by using nonscp- 
arable measurements there is a significant improvement 
in the accuracy of the estimation of unitary operations. 
It is also proven that in d dimensions the QCRB can be 
achieved if and only if the input state is maximally entan- 
gled. An open problem is the estimation of more general 
quantum operations, described by the Kraus decomposi- 
tion p|. 



putting things together 
d 2 -i 



and 



J2(k\t«) 2 = ^(<i>s\^ 



(23) 



(24) 
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